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Abstract 

It is widely believed that charm meson production and decay may play an impor- 
tant role in high energy astrophysical sources of neutrinos, especially those that are 
baryon-rich, providing an environment conducive to pp interactions. Using slow- 
jet supernovae (SJS) as an example of such a source, we study the detectability of 
high-energy neutrinos, paying particular attention to those produced from charmed- 
mesons. We highlight important distinguishing features in the ultra-high energy 
neutrino flux which would act as markers for the role of charm in the source. In 
particular, charm leads to significant event rates at higher energies, after the con- 
ventional (tt, K) neutrino fluxes fall off. We calculate event rates both for a nearby 
single source and for diffuse SJS fluxes for an IceCube-like detector. By comparing 
muon event rates for the conventional and prompt fluxes in different energy bins, 
we demonstrate the striking energy dependence in the rates induced by the presence 
of charm. We also show that it leads to an energy dependant flux ratio of shower 
to muon events, providing an additional important diagnostic tool for the presence 
of prompt neutrinos. Motivated by the infusion of high energy anti-electron neu- 
trinos into the flux by charm decay, we also study the detectability of the Glashow 
resonance due to these sources. 
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1 Introduction 



Over the coming decade, the detection of cosmic neutrinos is both an exciting prospect 
and a chahenge. The considerable efforts deployed towards its realization are made 
worthwhile by the potential to add to our knowledge, both of the astrophysics of high 
energy sources and of physics beyond the standard model [T]. The extant and forthcom- 
ing km^ sized facilities [21 [3l H] may detect high-energy neutrinos from sources such as 
Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB), heralding a new chap- 
ter in neutrino physics. Specifically, if such neutrinos are observed, their intensity and 
spectrum would provide key information towards understanding the dynamics of their 
source objects. They may also shed light on the origin of the observed cosmic rays. 
Additionally, flavor ratios for the three types of neutrinos may provide insights into new 
physics beyond the standard model, e.g., neutrino decay, the presence of pseudo-Dirac 
mass states and CPT violation [5l[6l[7]. The expected fiux ratio $e : : = 1 : 2 : 
in many sources is very sensitive to the existence of these effects. If a significant de- 
viation from the resulting democratic distribution $e : : '^'r = 1 : 1 ^ 1 expected 
at earth-based detectors (after oscillation effects are incorporated) [8] is reconstructed 
from the observations, it would be a clear signal of new physics. 

An unresolved issue which is expected to significantly affect the emitted and ob- 
served neutrino fiuxes and their ratios is the role played by charm meson production 
and subsequent decay. Such mesons are significantly short-lived compared to vr and 
K mesons, which are copiously produced in the jry and pp interactions which typify 
almost all such sources. This results in the charm meson neutrino spectrum in high 
energy astrophysical sources having qualitatively different features. Specifically, their 
decay is expected to lead to relatively flatter and elongated neutrino spectra compared 
to the sharply dropping contribution from vr and K decays. This leads to an enhance- 
ment of higher energy neutrinos in sources where charm plays a role, as shown recently 
in [9]. Ratios of fiuxes, which as stated above, are important diagnostics of new physics, 
would also be altered, since charm decay produces equal numbers of muon and electron 
neutrinos and antineutrinos. 

While these effects may not be important in all high energy sources, they could 
significantly alter the observed neutrino fluxes for a class of sources which are charac- 
terized by proton rich jets, or, in general, an environment which renders pp interactions 
important. It is our intent in this paper to examine the observational consequences for 
ultra high energy detectors for this class of sources, where charm mesons are produced 
in signiflcant numbers, along with vr and K. Our aim is to bring out general qualitative 
differences which would stand out as markers for the presence of charm production via 
pp interactions in the source. If observed, this would also shed light on the issue of the 
presence and relative importance of pp interactions versus p^ interactions. 

As a specific case in point, to facilitate quantitative estimates of event rates from 
which conclusions may be drawn, we choose as source slow-jet supernovae (SJS), which 
have been the focus of recent studies [9l [TOl [11]. However, several important features 
related to observations which emerge from our calculations for SJS may be generic to 
other types of sources whose jets are baryon rich. 

To this end, we calculate the event rates using the fiuxes given in [9l [ini [H]. This 
is done both for a single SJS source and for the diffuse flux obtained by summing over 
sources. While it was shown in [lOl [11] that a nearby source (at around 3 Mpc) is 
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observable via its neutrinos from vr and K decay, we show that the charm, or prompt 
component adds substantially to the neutrino event rate at energies significantly above 
those spanned by the conventional flux. Additionally, we find that due to the charmed- 
meson contribution, the diffuse flux from SJS alone can rise above the atmospheric 
neutrino background, and may already be constrained by the AMANDA observations [2]. 
It must be noted, however, that flux estimates of high energy sources, though usually 
underpinned by reliable physics, have inherent uncertainties. This is more so in the case 
of charm contributions, with their added QCD-related uncertainty. The event numbers 
we compute are intended to serve as estimates, whereas the qualitative observational 
features which are pointed out here, distinguishing the presence of charm contributions, 
may, in the long run, serve as more useful diagnostic tools as experiments accumulate 
data over the next decade. 

The plan of the paper is as follows: In Section 2, we briefly review the physics of SJS 
supplemented by charm meson production. We discuss the qualitative features of the 
flux and, for practical use, give its fitting functions. We then calculate the event rates 
and spectrum in an IceCube-like detector. In Section 3, we obtain the diffuse flux and 
associated event rate by summing up all contributions from SJS point sources. Section 
4 is devoted to a discussion of our results and conclusions. 

2 Neutrinos from Slow-Jet Supernovae 
2.1 The neutrino flux 

Since the discovery of the connections between long duration GRB and supernovae |12j . 
it has been conjectured that a significant fraction of core-collapse supernovae (SNe) 
may be accompanied by mildly relativistic jets which do not break through the stellar 
envelope. The choked jet may give rise to a radio after-glow without prompt 7 ray 
signals, and such objects may be detectable primarily only via their neutrino emission. 
The neutrino spectrum resulting from this scenario was first modeled by Razzaque, 
Meszaros, and Waxman (RMW) [10] and extended by Ando and Beacom (AB) [TT] . 
The essential difference between SJS and the related GRB jets is the mild gamma factor 
of the former, T ~ 3 and its baryon rich environment which makes pp collisions efficient. 

An extension of the work in [lOllllj has recently been discussed in [9]. These authors 
calculate the neutrino flux from charm meson production and decay in the source, by 
solving the cascade equations which describe the propagation of protons and mesons. 
The proton rich environment inside the slow jets is naturally conducive to charm pro- 
duction via pp interactions. Daughter neutrinos from charmed mesons typically have 
energies higher than those resulting from conventional ones. These mesons decay quickly, 
producing fg,^ -|- i^e,^i, before significant energy loss via interactions can occur, so that 
their contribution to the neutrino flux becomes dominant at > 10^'^ GeV. This result 
is a significant change in the overall features of the emergent neutrino flux and, to an 
extent, in its flavor composition. 

In general, meson spectra, whose shapes are inherited by the daughter neutrinos, 
are characterized by two power-law energy breaking scales. The first breaking point 
is determined by the competition between decay and hadronic cooling ( the latter 
characterized by, for instance, np and Kp interactions). The second energy IS 
determined by the transition between hadronic cooling and radiative cooling and marks 
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the onset of the latter. Below the first breaking energy, meson decay is dominant and 
the spectrum is the same as that of the shock accelerated protons, ~ E~'^. Between 
the first and the second scale, hadronic coohng is dominant, and the spectrum gets 
steepen to ^ E~^. Above the second breaking energy, radiative cooling dominates, and 
the spectrum falls as ~ E~'^ up to the cutoff given by the proton maximal energy. For 
pions, the two breaking energies are close to each other and relatively low; E^^^ ~ 30 GeV 
and -E^^^ ~ 100 GeV. Thus the pion flux retains an E~'^ power over TeV scales, up 
to the maximal cutoff energy. For kaons, the relevant breaking energies are given by 
~ 0.2 TeV and ~ 20 TeV [II]. 

For the charmed mesons and , due to their short lifetimes and large masses, the 
two scales are close to each other and have relatively high values, E^^^ ~ E^'^'^ ~ 10^ TeV 
[9]. Thus the charmed flux maintains a E'"^ power law all the way up to the proton 
maximal energy, leading to a significantly flatter spectrum. 

In Fig. [H we show various fluxes of muon neutrinos and antineutrinos at an Earth 
detector. The calculation uses RMW [lO], AB [llj and ERS [9] source fluxes and prop- 
agates them to the earth for an assumed SJS source at a distance of 20 Mpc. The fluxes 
are shown without neutrino oscillation effects, which we discuss in a later sub-section. 
We note that the flux estimates for muon neutrinos from vr and K differ in the three 
models, while exhibiting similar spectral shapes. The differences, in part, result from 
differing treatments of the maximal proton cutoff energies. 

For the ERS case, we also show the flux resulting from the incorporation of D^- 
meson flux (The curve for is similar to that for D^). We note that the ERS-D flux 
becomes dominant above Ei, ~ 10^ GeV. This flux, after propagation over a distance 
di from the source at redshift z is well represented by the following fitting function. 



Fd 



27r4 



(1 + ^) 



exp 



(1 + z)E, 



(2.1) 



where 



E, 



D 
br 



2.6 X 10^° GeV • s~ 
lO'^-^ GeV, 
1.35. 



(2.2) 
(2.3) 
(2.4) 



Similarly, for the conventional ERS [tt + K meson) flux, we find 
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where 



Egr 



Pc 



3.1 X 10^° GeV • s" 




Eu < Eg, 



br' 

for E,>Eg. 



(2.6) 
(2.7) 

(2.8) 



These expressions are useful not only for calculating the event numbers but also to 
deduce the diffuse flux, both of which we discuss below. The above fluxes do not 
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accommodate the effect of the meson accelerations by the internal shock [T3]. We note, 
however, that accounting for this acceleration and subsequent losses due to the increased 
lifetime and increased interactions would lower the ERS flux estimates, and the event 
rate predictions we make in this paper for them. 




2.2 Neutrino oscillation effects 

The neutrinos produced at a source can change flavors during propagation to earth. 
Although their energies are very high, the distance between a source and a detector is 
also large enough to average out the oscillatory behavior. Writing the vacuum transition 
probabilities between flavor eigenstates a and (3 in terms of the lepton mixing matrix V 
only, we have 

-fa/3 = iKfiPlVaiP 

s s \ 

Ul-s) i(l-s) , (2.9) 
\{l-s) \{l-s)l 

where s = cos^ 9i2 svc? 9i2 and 612 is the solar angle in the standard parameterization for 
V. The simple form of Eq. (j2.9p is a result of setting the "reactor" and "atmospheric" 
mixing angles as ^13 = 0° and 623 = 45° respectively. As a consequence of this maximal 
fj,-T mixing, the flux ratio of and at earth must be 1. In addition, the standard 
ratio $e '■ : $r = 1 : 2 : at a source means results in a democratic configuration 
$e : : = 1 : 1 : 1 at a detector, independent of the value of 612 ■ 

For SJS, the environment is baryon rich and muons lose their energy before decaying. 
Thus, for the pion and kaon fluxes, the initial configurations are $e '■ : $r = : 1 : 0. 
On the other hand, charmed mesons decay semi-leptonically such that the initial ratio 
is $e : : $T- = 1 : 1 : 0. According to Eq. (|2.9p . these initial fluxes at the source 
win be transmuted to P ■ (0, 1,0)'^ = (s, (1 - s)/2, (1 - s)/2f = (0.22, 0.39, 0.39)^ and 
P- (1,1,0)'^ = (1 - s, (1 + s)/2, (1 + s)/2)T = (0.78,0.61,0.61)'^ at earth, where we 
have used 612 = 34.5° as a current best fit value [H]. To summarize, we note the fluxes 
Eq. (|2.ip and Eq. (|2.5p after oscihation are given by = s • Fc, = (1 — s) • Fd, 
F^ = (1 — s)/2 • Fc, and F^ = (1 + s)/2 • Fd, denoting the flavor species by upper 
subscripts. 

The difference of the initial flavor ratio between the prompt and the conventional 
fluxes and their different spectral behavior gives rise to an energy dependent flavor ratio 
at earth. This is given by 

^ = 2 + ,2 10) 

where ^c,D stand for generic charm and conventional fluxes. The ratio is bounded as 
follows, 2s/(l -s) K. 0.56 < ^e/^fi < 2(l-s)/(l + s) 1.28 by the high and low-energy 
asymptotic behavior, marked by the dominance of the prompt and the conventional 
component respectively. The energy dependence is most significant around the energy 
at which the prompt and conventional components are equal. We show the behavior of 
the ratio Eq. ()2.10p by the solid curve in Fig. [2l by plotting the prompt and conventional 
flux Eq. (j2.ip and Eq. (j2.5p at (i^ = 20 Mpc. Of course, this remains unaltered when we 
compute the diffuse flux by integrating over the distribution and luminosity of sources. 
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Table 1: The number of up-going ly^N — > fj,X events from an SJS at 3 Mpc. The 
columns show the events in each decade of energy (in GeV) and the total events. An 
instrumented volume of 1^ = l(km^), and an energy threshold i?™™ = 10^ (GeV) is 
assumed. 

since it depends only on the mixing matrix and the fact that the oscillation length is 
much shorter than the source distance. 

Fig. [2] highlights the importance of detecting both electron and muon flavored neu- 
trinos, particularly in the high-energy region above 50 TeV. The rise of the ratio above 
unity at high-energy will be an indication that charmed mesons contribute to neutrino 
production. In addition, the observation of the transition of the ratio around 1 ~ 100 
TeV would indicate that the source environments are baryon rich, to a degree that allows 
the charm component to rise above the conventional tt, K component and the muon to 
be damped. We also note the significant difference in the SJS without charm and SJS 
with charm curves in the figure, i.e from a constant 0.5 for the former, to an energy 
dependant transition from 0.5 to 1.25 for the latter. Thus, shower events would register 
an increase from their low values compared to muon events at low energies, and rise 
above them at high energies. 

We note that the canonical flavor ratio <I>e : <l>;i : $t = 1 : 2 : 0, the ratio for 
muon damped sources, <l>e : <I>^ : <5t = : 1 : and that for semi-leptonic D decay, 

■ ■ = ^ ■ ^ '■ should all be regarded as idealizations which can be used as first 
approximations. More precisely, there can be small deviations from the above standard 
ratios with some uncertainties [15]. For example, it has been discussed that the muon 
damping cannot be perfect and z^e component consequently never can fall to zero |15j . 
Inclusion of these effects will amount to finite thickness for the curves presented in 
Fig. [2j However, qualitative behavior such as the transition of the ratio around 1 ~ 100 
TeV will not change significantly and the curves in Fig. [2] still represent the essential 
qualitative physics of the charm effect. 

2.3 Detection of high-energy neutrinos 

We next compute the up and down-going event rates for an IceCube-like detector for 
the fluxes considered in the previous sub-section, assuming a point source at a distance 
of 3 Mpc. While this has been previously done in [10\ [TT] for the conventional flux, 
our purpose here is to compare this contribution quantitatively with that of the charm 
mesons. 

For up-going events, we incorporate the attenuation of neutrino flux by the earth 
matter, and also consider partially contained events which enhance the effective volume 
of the detector significantly for muon events. The event rate is given by [16] 

Rate = ANa fdE, ■ {RiEf^^, E,)) ■ acc{E,) ■ S{E,) ■ F^ ,, (2.11) 
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Table 2: Same as Table 1, but for down-going i^^N — s- fiX events. 

where A is the effective area of the detector, and is Avogadro's number A^^ = 
6.022 X 10^'^ (cm~^). The muon range {R{E^™,Ei,)) stands for the average distance 
traversed by a muon in earth matter. Muons created by neutrinos with energy Ey 
attenuate to the energy E"™™, after traveling {R{E™™^ E^)). The factor S{Ey) represents 
the shadowing (i.e attenuation) of the up-going neutrinos by the earth. 

For down-going events, the event rate for the vN — > jiX process is given by 

Rate = NAVes J dE^ ■ acciK) ■ F^ j^ (2.12) 

where Vcs is the effective volume of the detector, assumed here to be 1 (km^). cjcc is 
the neutrino-nucleon charged current cross section. 

In addition to the vN fiX events, we have also calculated events for the Glashow 
resonance process, i^eC — > something, which peaks around E^, = 6.3 PeV, primarily 
because we expect an enhanced contribution from charm at the higher energies. The 
event rate for down-moving neutrinos from this channel is given by 

10 f 1 

Rate = —NaV,s I dE, ■ ap^e{E,) ■ -F^^^ (2-13) 

Our results for all the various event spectra are shown in Fig.[3ljl|) for up(down)-going 
events. Clearly, the charm contribution dominates above ^ 50 TeV, and its contribution 
is significant even around 500 — 1000 TeV, where the vr and K contributions to the event 
rate are negligible. Table 1(2) shows the integrated number of events for each energy 
decade with up(down)-going muon events. First, we note that the ERS flux, even for 
the conventional case, yields rates which are about 20-30 times higher than the AB or 
RMW fluxes, due to an enhanced pp contribution in this model. The charm contribution 
for ERS is a full third of the integrated rate for the conventional flux. Importantly, it 
manifests itself in energy regions where the conventional contribution is low, i.e., in 
the band lO'^ — 10^ GeV. This both enhances its prospects for detection and acts as a 
distinctive signature for charm. Finally, we note that since the typical duration of the 
burst is a few seconds, the background from the atmospheric flux is small. 

In Fig. m we show the event spectrum for i>f,e — > v^^i channel as well. From the 
figure, we can clearly see the resonant behavior at 6.3 PeV. The number of events for 
this mode is, with the ERS-D flux, 0.06 in the 10^~^ GeV bin. While this contribution is 
~ 50 % of its total event tally over the full energy range of 10^~® GeV, it is unfortunately 
too tiny to be detectable. 
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3 The Diffuse neutrino flux from SJS 



3.1 The diffuse neutrino flux 



Prior to presenting our results for the events from the diffuse SJS fluxes, we note that 
charm is likely to contribute to the atmospheric flux as well. This has been discussed 
in detail in, for example, in [T7|. We have ascertained that its contribution is low and 
does not constitute a significant background to the SJS flux. 

The diffuse neutrino flux is obtained by integrating the point source flux over all 
SNe, weighting the flux by a rate the number of SNe events per unit time per 



solid angle covering the earth sky. The distribution ^j^^ is given by [10 



dtdn 

where the SNe rate hsn{z) is 

nsn(^) = 0.017 



(1 + ^) 



dt 
dz 



0.32e^-^^ 



,3.82 



+ 45 



X 10" 



-81 



cm 



The cosmic time t and z is related as 

dt 1 1 

dz 



Hol + Z v/(l + nmZ){l + Z)2 - nA{2z + 



(3.1) 



(3.2) 



(3.3) 



where = 0.3 and f^A = 0.7 within the standard cosmology. With these ingredients 
at hand, the diffuse neutrino flux is given by (e.g. for charm flux) 



^diff 



dz ■ 
dtdi I 



4^2 R2 



dz ■ , 
(1 + 



dt 






dz 


exp 


;( 



'{\^z)E^ 



(3.4) 



where tj is the jet duration of ~ 10 s, and represents a rate of SNe with an accom- 
panying jet. The probability that the jet points to the earth is assumed to be l/2r2. 

In Fig. O we show + diffuse fluxes. From the figure, we can see that the 
total contribution for the ERS flux already begins to breach AMANDA data, and hence 
the charm contribution is likely to be both tested and constrained over the coming 
years. As already pointed out, the range over which it is dominant is largely different 
from that spanned by the conventional contribution from vr and K mesons, making it 
observationally distinct. 

Table 3 shows the number of the events per year per steradian for each flux. In Fig. [S] 
and Table 3, we have, somewhat unrealistically, endowed all supernovae with slow jets 
of r = 3, and thus taken = 1. For the ERS fluxes, the diffuse rates are quite high 
and certainly observable, especially in the higher background free energy bins. Roughly, 
a total of 200 up-going and 180 down-going events are expected per yeaJl. A more 
realistic value of ^sn would lower our rates by an appropriate factor. 



We note that the Waxman-Bahcall bound [18], which appears to be breached, does not apply to 
such sources since they do not contribute to the cosmic-ray flux 
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Table 3: The number of up-going muon events per year per steradian with the diffuse 
prompt and conventional fluxes. The instrumented volume of the detector is taken 
as y = 1 (km^); the effective area A = 1 (km^) for up-going and the effective volume 
= 1 (km^) for down-going. We are taking the energy threshold as E"^™ = 10^ (GeV). 
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1 














1 


RMW-K, TT 
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Table 4: Same as Table 3, but for down-going muon events. 



Finally, we present simple fits to the flux curves for the muon neutrinos shown in 
Fig. [5l The fluxes are well represented by 

y = ax'' + ce'^'' , (3.5) 

where y = logiQ{^El),x = logiQ{Eu). The diffuse flux <I> and the neutrino energy E^ 
are in the same units as in Fig. [5j Using Eq. (j3.5p . the neutrino flux for each case, (i) 
SJS conventional (ii) SJS L>-meson and (iii) SJS total can be fitted and the sets of the 
best-fit parameters, a, b, c, and d are listed in Table El 

4 Summary and Conclusion 

The production and decay of charmed mesons can play an important role for the de- 
tection of high-energy cosmic neutrinos. Slow-jet supernovae provide both an excellent 
example and a test case in which the charm effect can be significant due to the baryon 
rich environment, leading to the rise oi pp contributions to the total neutrino rate. This 
enhances the detectability of the SJS neutrinos not only with the point sources but also 
with the diffuse fiux. 

The charm contribution appears to open up the possibility of detecting SJS via their 
diffuse flux, using the ERS flux calculations. While the absolute intensity of the fiux 
should be regarded with caution at the present stage, the elongated shape of the total 
flux offers a shot at detection, provided pp interactions play a significant role, and if the 
the jet parameters (and the faction ^sn ) are favorable. 

By virtue of the semileptonic decay of the charmed mesons, the flavor ratio $e/^/i 
must be ~ 2(1 — sin^0i2), which is essentially different from the typical AGN or GRB 
case. Furthermore, the energy dependence of ^^/^fj^ also carries potentially important 
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a 


b 


c 


d 


(i) SJS conventional 


-0.86203 


1.41904 


-25.73761 


-0.84254 


(ii) SJS D-meson 


-6.20429 


0.06759 


-6.8552 X 10-y 


2.83302 


(iii) SJS total 


-5.17119 


0.16516 


-2.1531 X 10"^ 


2.68061 



Table 5: Best-fitted values of the parameters a, b, c, d for different SJS fluxes. 



information. The ratio shifts occurs around the crossing point i^cross where the D- 
meson flux becomes equal to the kaon flux. The crossing energy is written as -Ecross ~ 
(/x//d)^x^ « 200^i^\ where /k(/d) ^ 0.06(3 x lO""^) are the kaon and the charmed- 
meson multiplicities, and is the first breaking energy for kaon. The breaking energy 

E^^^ is very sensitive to the bulk Lorentz factor F such that E^^ oc r^*^^ depending on 
the assumption on the jet-opening angle [TT] . 

To summarize, charmed mesons can play an important role in a certain class of 
high energy sources in general and some supernovae in particular, where jets are mildly 
relativistic and environments baryon-rich. Detection of the charm related neutrino flux 
may be possible with current and near- future detector sensitivities, and can offer a 
valuable handle on several important open questions related to the dynamics of such 
sources. 
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Figure 1: Muon neutrino and antineutrino fluxes from a source at 20 Mpc. The 
thick(thin)-solid curve shows the D^{K^ + vr^) flux in [9]. The dashed and dotted 
curves are the + vr^ fluxes in [10] and [H] respectively. The fluxes are given in the 
earth-observer frame. The effects of neutrino oscillation are not incorporated in this set 
of curves. 
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Figure 2: The e/ fx flavor ratio of the flux, as a function of the neutrino energy. The solid 
curve shows the full contribution with SJS. The dotted line is for SJS, but only with 
the conventional component. The dashed hue represents a typical value for the sources 
which allow full muon decay, such as AGN (without charm). 
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Figure 3: Event spectrum of up-going i/N fi~X with a source at c?i = 3 Mpc for 
ERS flux [9]. The blue (thick solid), red (thin solid) and the dashed curve shows D^, 
K and tt contribution respectively. 
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Figure 4: Same as Fig. [3] but for down-going events. The lower curve in each case is for 
z^eC — > z^^/i channel. 
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Figure 5: The diffuse i^^ + z^^j flux from slow-jet supernovae, compared with the recon- 
structed atmospheric neutrino flux by AMANDA 2000 year data [2]. The thick plain 
(dashed) curve shows the total (conventional) flux for SJS, and the thin dashed curve is 
the prompt flux from the atmosphere [TTj. The horizontal dashed lines shows the WB 
limit ~ 1 - 5 X 10^8 GeV • cm-^ . g-i • sr-^ [l8]. 
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